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The morphology of galaxy clusters is quantified using Minkowski functionals, espe- 
cially the vector— valued ones, which contain directional information and are related 
to curvature centroids. The asymmetry of clusters and the amount of their sub- 
structure can be characterized in a unique way using these measures. — We briefly 
introduce vector-valued Minkowski functionals (also known as QuermaB vectors) 
and suggest their application to cluster data in terms of a morphological charac- 
terization of excursion sets. Furthermore, we develop robust structure functions 
which describe the dynamical state of a cluster and study the evolution of clusters 
using numerical simulations. 



1 Introduction 

The substructure of galaxy clusters is important not only for its own sake. For 
the usual mass estimates of clusters depend on how the mass distribution is 
modelled, especially whether substucture is ignored or not. Furthermore, the 
amount of substructure observed in clusters can serve as a sensitive probe of 
the cosmological parameters, especially the overall matter density ^o^- 
In this context it is necessary to measure cluster substructure in a quantitative 
way. So far, only very crude substructure measures have been developed (for 
an overview see Crone et al.B, for a different approach Grebenev et al.u). In our 
contribution, we use Minkowski functionals to characterize cluster morphology 
in the systematic framework of integral geometry. The Minkowski functionals 
can be introduced in an axiomatic way (Section-M) and have been sucessfully 
applied in discriminating large scale structure oil We focus on vector-valued 
Minkowski functionals which also incorporate directional information. 



2 The Minkowski functionals 

The (scalar) Minkowski functionals are defined for closed bodies and can be 
characterized by simple requirements such as motion invariance, addivity and 
conditional continuity (see Figure ^) . Although these requirements are fairly 
general, Hadwiger's theoremQ states that in d dimensions there are only d + 1 
linear independent Minkowski functionals Vi . This means that a characteriza- 
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Motion invariance: Additivity: Conditional continuity: 

V(A) = V(gA) V(AuB)=V(A)+V(B) - V(A nB) V(A n )->V(A) 



Figure 1: The requirements defining the Minkowski functionals 

tion of bodies is possible in a unique way in terms of d + 1 measures which 
carry intuitive meanings. For example, in two dimensions, the first Minkowski 
functional, called Vq, is simply the surface content. The other two Minkowski 
functionals are the length of the circumference and the Euler characteristic; 
they can be obtained by integrating over the Lebesgues surface of the body 
and weighting with curvature measures. The Eulfir characteristic V% is related 
to the genus, for further details see Mecke et al.B. 

An interesting extension can be made by going to vector-valued Minkowski 
functionals. These transform like vectors (motion egrovariance), and it turns 
out that to each scalar Minkowski funtional Vi there corresponds a vector- 
valued brother Vi or equivalently a normalized (curvature) centroid Pi = -pf, 
which is simply a point within .-the convex covering of the body, for further 
details compare Beisbart et al. Q. In Figure || some examples are given how 
Quermafi vectors work. Most notably, the (curvature) centroids are sensitive 
to symmetry. For a roughly symmetric body the centroids coincide, while they 
split up for an asymmetric configuration. 

3 The method: Minkowski functionals of excursions sets 

In order to employ the Minkowski functionals for the quantification of cluster 
morphology we have to associate bodies with the usual data sets (2d galaxy/X- 
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Figure 2: Examples 
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Figure 3: The Minkowski functionals Vq, V±, V2 and the morphological parameter £ are 
shown as functions of the density threshold u. 

ray-photon positions). For this purpose we smooth the galaxy positions with 
a Gaussian kernel to obtain a density field u(x). The smoothing length serves 
as a diagnostic parameter determining at which scale substructure is resolved. 
The morphology of the excursion sets M v = {x s.R d | u(x) > v\ can now be 
described in terms of the Minkowski functionals a For spherically symmetric 
reference models the excursions sets become simply circles, and the curvature 
centroids Pi coincide. 

This changes if we proceed to analyzing more realistic data. For example, we 
investigate clusters simulated by Matthias Steinmetz for an f2rj = 1-universetJ. 
In Figure || the behaviour of the scalar Minkowski functionals Vi is shown in 
dependence of the density threshold defining the excursion set; so we get in- 
dividual signatures of the cluster: the Euler characteristic V2 simply counts 
the components of the cluster; for high density thresholds only the core of the 
cluster remains (resulting in V% = 1), whereas for lower levels up to 10 com- 
ponents are counted; the second Minkowski functional V\ (second panel), the 
length of the circumference, gives an averaged, smoothed and inversed density 
profile of the cluster. The ratio of the first two Minkowski functionals may 
serve as a measure of how crooked the isodensity curves (the surroundings of 
the excursions sets) are. - The wandering of the centroids withj-the density 
threshold is a natural extension of the well-known centroid shift □. However, 
to quantify the amount of substructure present in the cluster it is better to 
compute the surface content C, and the length of the circumference ij of the 
triangle formed by the (curvature) centroids. The behaviour of the so defined 
morphological parameter £ can be seen in the fourth panel of Figure |j| 
To further condense the detailed information present in the Minkowski func- 
tionals we form parameters a, (3 and 7, which overcome the dependence on 
the density 0. For example we define a by integrating the deviation of the 
Euler characteristic from 1 , which is the value for a highly symmetric reference 
model. All these parameters are constructed in such a way that they assume 
the value if the cluster is spherically symmetric. 
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Figure 4: Evolution of simulated clusters (Qp = 1) m the phase diagram 



4 Results and further prospects 

Using the above defined morphological descriptors we can follow the dynami- 
cal evolution of simulated □ clusters within a phase-diagram defined in the 
a/?7-space. Two examples are shown in Figure ^|. The paths of evolution take 
similiar tracks towards a more relaxed, substructure-poor state. Note, however, 
that these tracks span different scales: cg02 has a higher amount of substruc- 
ture than cg04. In different cosmogonies the paths will differ; e.g. in low-fio 
universes the morphological evolution will stay almost within the /37-plane. 
Hence, the cosmological parameters can be determined quantitatively in com- 
parison with observed cluster data. Projection effects in the substructure may 
be controlled using the stereological properties of the Minkowski functionals. 
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